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Abstract 

Behavior of planar multilayer periodic structures due to plane wave excitation has been 
studied using the transfer matrix method. Multilayer structure is taken with periodicity 
two. That is, layers at even and odd locations repeat themselves. Layers at odd locations 
are of chiral nihility metamaterial whereas three different cases for layers at even locations 
are considered, i.e., dielectric, chiral and chiral nihility. Effects of polarization rotation 
due to the optical activity is studied with respect to the angle of incidence and frequency 
in terahertz domain. Chiral nihility introduces property of transparency to the structure 
for normal incidence while complete rejection is observed for chiral nihility-chiral nihility 
structure at oblique incidence . 
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0.1 Introduction 



Objects in our universe can be classified in two groups: achiral and chiral. An object 
not superimposable to its mirror image by any translation or rotation is called chiral ob- 
ject. Achiral object are non-chiral and can be superposed to its mirror image [1, 2]. A 
metamaterial microscopically composed of chiral objects is called chiral metamaterials. The 
constitutive relations for isotopic, homogenous and non-diffusive chiral material are written 
below [3]: 

D = eE+(x-i/^)^/^H (1) 
B = /iH + (x + jK)V^E (2) 

where Hi and x chirality and non-reciprocity parameter respectively. Moreover e and /i 
are permittivity and permeability representing the electric and magnetic polarizability of 
the material. 



Lakhtakia proposed the concept of nihility metamaterial and nihility means real parts 
of permittivity and permeability are zero [4] . He studied the scattering of electromagnetic 
plane waves from a cylinder composed of nihility metamaterial [5]. Later, Tretyakov et 
al. extended the concept of nihility for the isotropic chiral metamaterials [6]. Chiral 
nihility (CN) is a special case of chiral metamaterial for which real parts of permittivity 
and permeability are simultaneously zero at certain frequency, i.e., e — > and /i — > 0. 
The constitutive relations for CN metamaterials are [6-8] , 

D = -jf^^yeofioH (3) 
B = j/^y^E. (4) 

Several researchers have studied the behavior of metamaterial interfaces and structures 
when these are exposed to electromagnetic excitation. Some special arrangements using 
metamaterials had been proposed to fabricate absorbers, antenna radomes and cloaks [9- 
11]. Interaction of electromagnetic waves with chiral slabs with planner interfaces had been 
studied for different apphcations [12-14]. Rejection and tunnehng of electromagnetic fields 
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from chiral and/or CN planar interfaces and waveguide composed of planar CN metamate- 
rials had also been studied [15, 16]. Reflection and transmission from CN slabs embedded 
in other materials or backed by fractional dual interfaces are available in [17, 18]. Khalid 
et al. had studied scattering of electromagnetic waves from cylindrical DB boundaries in 
the presence of chiral and CN metamaterials [19]. Taj et al. had studied the chiral and 
CN metamaterials interfaces as focusing surfaces by using Maslov's method [20] . Sabah et 
al. had studied the behavior of multilayered structure. Sabah along with co-workers had 
calculated the response of the chiral mirror by using transfer matrix method (TMM) [21]. 
Polarization rotator by using multilayered structures had been studied by using conven- 
tional chiral and chiral metamaterial for normal incidence [22] . Planner interfaces composed 
of chiral metamaterials had studied for oblique incidence as filters in [23, 25]. 

In current paper, planar periodic multilayered structures having odd number of slabs 
are analyzed. Two layers form a period and it is assumed that layers at odd locations are 
of CN metamaterial. Three different cases are considered for layers at even locations, i.e., 
dielectric, chiral, and CN. Purpose of the study is to find the effects of CN metamaterial 
on refiected and transmitted powers. 



. 2 Formulation 



Consider planar periodic stratified structure of chiral metamaterials as shown in Figure 1. 
Each slab is of infinite extent in x and y directions of the Cartesian coordinate system. 
Slabs at odd locations have identical constitutive parameters while slabs at even locations 
have same values of constitutive parameters. Therefore, each slab at odd location is labeled 
as A whereas each slab at even location is labeled as B. It is also assumed that slab labeled 
A have higher value of refractive index than slab labeled B. Constitutive parameters for 
slab A are denoted by (eij,/iij, kh) and constitutive parameters for slab B are denoted by 
(e^, /iL, i^l)' The stratified structure is placed in air medium having constitutive parameters 
(eo, /io). 
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Figure 1: Stratified metamaterials at oblique incidence 

The structure is excited by a parallel polarized monochromatic plane wave. The expres- 
sion for incident electric field is given below: 

E^ = [Ei\\{xcos9i + zsin^^)] exp[— j''/co(^cos^^ — xsin^^)] (5) 

where /cq = Usje^ji^ and 9i is the angle of incidence with respect to z-axis. 

It is assumed that refiected and transmitted fields contain both parallel and perpendic- 
ular components. Expressions for refiected and transmitted waves, in terms of unknown 
coefiicients, are written as 

E^ = [£'^||(xcos^^ — zsin^^) + Eri_y\ exp[— j''/co(— ^cos^^ — xsin^^)] (6) 
Et = [Et\\{xco^6t — z^iiiOt) + Et±_y\ exp[— j''/co(^cos^t — xsin^t)] (7) 

Or and 6t are angles of refiection and transmission in air with respect to z-axis. 

Electric and magnetic fields inside the m-th chiral slab are written as linear combination 
of left circularly and right circularly polarized plane waves. Both waves are propagating 
in forward as well as backward directions, that is, in positive and negative z-directions. 
Expressions for total electric field in m-th chiral slab are given below [9]: 

E+ = E+^exp[-j/c^L(^cos6>^L - xsin6^^L)] + E+^exp[-j/c^i^(z cos (9^i^ - xsin6>^i^)], 
E^ = E;;^exp[-j/c^L(-^cos(9^L - xsin6>^L)] + E;;^exp[-j/c^i^(-z cos 6^^i^ - xsin6^^i^)]. 
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with 



^mL = E+^{xCOSemL + Z sin e^L+jy), 

^mR = E^r{^ cos OmR + Z sill OmR " jv) , 

'^mL = -^mL sin 6'^^ - X COS 9^1 +jy), 

^mR = ^mRi^ ^ill 0^r - X COS 6,riR " j^) • 



(8) 
(9) 
(10) 

(11) 



Superscripts + and - are used to describe the waves propagating in forward and backward 
directions, respectively. Moreover Oj^l and 9jnR are used to describe the angles of LCP and 
RCP waves with respect to the z-axis. The magnetic fields outside and inside the chiral 
slab may be obtained from Maxwell curl equations. kmR and kmL are wave numbers for 
LCP and RCP waves, respectively and are given as. 



£"^11, Et\\^ Er±^ Et±^ E^j^ and E^j^ are unknown coefficients to be determined using the 
boundary conditions. Transfer matrix method (TMM) is used to determine the reffected 
and transmitted fields from the planar stratified structure. First of all, matching matrices 
are calculated by imposing the boundary conditions at each interface. To relate fields at two 
consecutive interfaces, propagation matrices are also computed. Product of these matrices 
finally relates the reflected and transmitted flelds with the incident fleld in form of matrix 
equation as given below [25], 



(12) 
(13) 



kmR = u;{i^i^/eofio + yf^i), 




(14) 



where relating transfer matrix is given below, 



Mi.Pa.T^.M2 
Mab-Pb-Mba.Pa 



(15) 
(16) 
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with Ml and M2 matching matrices which relate fields across first and last interface of 
planar stratified structure, respectively. In general, matching matrices Mij correspond to 
an interface with slab i on left side and slab j on right side of the interface. Propagation 
matrices are represented by Pa and Pg corresponding to slab A and slab B, respectively. 



0.3 Results and Discussions 

In this section, numerical results for three different structures are obtained to note the 
effects of CN metamaterials on transmitted and refiected powers. For all cases, parallel 
polarized plane wave {Ei_\_ — 0) is considered for incidence and structure have odd numbers 
of slabs. In all cases, each slab have optical width Ao/4 except for CN case when slabs 
have physical width equal to Ao/4, where Aq corresponds to the first harmonic frequency 
of /o =lTHz. 



0.3.1 CN-dielectric Structure 

In the first case, slabs labeled as A are of CN metamaterial whereas slabs labeled as B 
are of dielectric medium. It is assumed that structure is composed of five slabs. Each CN 
slab has physical width of Ao/4 and each dielectric slab have thickness Ao/4n^, where rid is 
refractive index of the dielectric slab. Figure 2 describes behavior of refiected power versus 
frequency for normal incidence. It has sinusoidal behavior with maximum amplitude equal 
to 40 percent of the incident power. Refiected power approaches to zero at 1.0, 2.0, 3.0 and 
4.0THz frequencies. Figure 3 shows the behavior of transmitted power versus frequency. 
Transmitted power has both cross and co-polarized components and these are shown by 
dotted and solid lines, respectively. It is also observed that rotation is 45^ at 1.0 and 
3.0THz, and complete rotation is noted at 2THz for normal incidence. 

In Figure 4 and Figure 5, refiected and transmitted powers versus frequency for oblique 
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incidence are plotted. It is noted that total reflection is observed at higher frequencies 
but some ripples are also noted in the reflected and transmitted powers at lower values of 
frequencies. 

Figure 6 shows reflected power versus angle of incidence. Total reflection is noted for 
6i > 22^ for flxed frequency of ITHz. The increasing trend for co-polarized component is 
also observed and it takes value unity from 6i — 60^ to 9i = 90^. Corresponding transmitted 
power, shown in Figure 7, has both cross and co-polarized components up to 6i — 22^ having 
maximum rotation at 15^ and becomes zero for next values. 

0.3.2 CN-CN structure 

For second case, it is assumed that both slabs A and B are of CN metamaterials with 
same value of chirality but having different limiting values of permittivity and permeabil- 
ity. Physical width of each slab is equal to Ao/4. In Figure 8 and Figure 9, it is observed 
that total power is transmitted through the structure but due to optical activity of the 
chiral material, the plane of polarization gets rotated. Rotation is 90^ at 2THz frequency 
for the chirality parameter hch — k.l—O.I 

In Figure 10 and Figure 11, behavior of CN-CN structure is studied versus frequency 
for oblique incidence. For these flgures, 9i — 45^ is considered and total reflection is noted 
after fraction of ITHz and total transmission is noted for 0.2THz. 

Next two flgures describe behavior of reflected and transmitted power versus angle of 
incidence. In Figure 12, complete reflected power is noted from 9i — 15^ to 9i — 90^ having 
maximum contribution of co-polarized component and cross component has largest contri- 
bution at 9i = 10^. In Figure 13, transmitted power is observed only at lower values of 
angle of incidence with major contribution from co-polarized component. 
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0.3.3 CN-chiral structure 



In third case, we consider that the slab at odd positions are of chiral metamaterial. Chiral 
slabs have optical width equal to quarter wavelength and lower chirality value than CN 
slabs. In Figure 14 and Figure 15, it is noted that reflected power has small percentage 
of incident power and transmitted power has major contribution of the cross polarized 
component for small values of angle of incidence. Whereas at higher degrees of incidence 
angle reflected power starts increasing and transmitted power shows decaying trend. 
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Figure 2: en = 1-6 * 10"^ = 1 * 10-^ ul = 2.2, cIh = Inildi = Ao/4, kh = 0.167 
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Figure 3: ch = 1.6 * 10-^ fin = I ^ 10"^ Ul = 2.2, rf^^ = \nL\dL = Ao/4, /^^^ = 0.167 
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Figure 4: e^f 1.6 x 10 ^, /x/f = 1 x 10 ^, 2.2, rf^f I^IlI^^l Ao/4, hiH = i^l = 



0.1, = 70^ 
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Figure 5: en 1.6 x 10 ^, /x/f = 1 x 10 ^, tll 2.2, dn Ao/4, = i^l = 

0.1, = 70^ 
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Figure 6: ch = 1.6 x 10 ^, /xj^ = 1 x 10 ^, Ul = 2.2, dn = I^IlI^l = Ao/4, k.h = 0.1, 
f/fo = 1 
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Figure 7: ch = 1.6 x 10"^, /x^f = 1 x 10"^ Ul = 2.2, rf^^ = \nL\dL = Ao/4, /^^^ 
///o = 1 
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Figure 8: en = 1.6 x 10-^ = 2.5 x 10-^ dn = di = Ao/4, kh = kl = 0.1 
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Figure 9: ch = 1.6 x 10-^ cl = 2.5 x 10-^ dn = dL = Ao/4, i^h = i^l = 0.1 
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Figure 10: en = 1-6 x 10-^ = 2.5 x 10-^ = = 10-^ dn = di = Ao/4, 
kh = kl = 0.1, = 45 
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Figure 11: ch = 1.6 x 10-^ cl = 2.5 x 10-^ fiH = I^l = 10"^ dn = dL = Ao/4, 
i^H = f^L = 0.1, 9i = 15 
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Figure 12: en = 1-6 x 10-^ = 2.5 x 10-^ = I^l = 10-^ dn = di = Ao/4, 
I^H = hl = 0.1, ///o = 1 
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Figure 13: en = 1.6 x 10-^ cl = 2.5 x 10-^ dn = (11 = Ao/4, fin = I^l = 10"^ 
i^H = i^L = 0.1, ///o = 1 
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Figure 14: ej^ = 1.6 x 10-^ = 2.5 x 10-^ i_ih = I^l = 10-^ dn = cLl = Ao/4, 
kh = K,L = 0.1, ///o = 1 
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Figure 15: ch = 1.6 x 10-^ cl = 2.5 x 10-^ dn = dL = Ao/4, fiR = I^l = 10"^ 
i^H = i^L = 0.1, ///o = 1 

0.4 Conclusions 

Behavior of planar multilayer periodic chiral nihility structures due to plane wave excitation 
has been studied using the transfer matrix method. Each multilayer structure is taken with 
periodicity two. First of all, behavior of CN-dielectric structure has been studied for the 
case of normal incidence. Total transmission as well as maximum rotation of power happens 
at 2THz. This stratified structure shows pronounced effects for oblique incidence. Total 
refiection of power for higher frequencies and small ripples, both in refiected and transmitted 
powers, for lower frequencies at oblique incidence are noted. Refiected power approaches 
to unity by increasing the value of angle of incidence and corresponding transmitted power 
has major part for small values of incidence angle. 

Complete transparent behavior of CN-CN structure for all frequencies has observed for 
normal incidence. Chirality introduces polarization rotation having maximum rotation 
at even harmonics. But for oblique incidence, total refiection from same structure has 
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been observed. At the end, behavior of CN-chiral structure has been studied and total 
rejection is observed for small incidence angles having major contribution of cross polarized 
components. 
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